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Abstract 


A batch air combat simulation environment , the tactical maneuvering 
simulator (TMS), is a tool for developing and evaluating tactical ma- 
neuvering logics that can also be used to evaluate the tactical implica- 
tions of perturbations to aircraft performance or supporting systems. The 
TMS can simulate air combat between any number of engagement par- 
ticipants , with practical limits imposed by computer memory and process- 
ing power. Aircraft are modeled using equations of motion , control laws , 
aerodynamics , and propulsive characteristics equivalent to those used in 
high-fidelity piloted simulation. Data bases representative of a modern 
high-performance aircraft with and without thrust- vectoring capability are 
included. To simplify the task of developing and implementing maneuver- 
ing logics in the TMS , an outer-loop control system, the tactical autopilot 
(TA), is implemented in the aircraft simulation model. The TA converts 
guidance commands by computerized maneuvering logics from desired an- 
gle of attack and wind-axis bank angle to inputs for the inner-loop control 
augmentation system of the aircraft. This report describes the capabilities 
and operation of the TMS and the TA. 


Introduction 

As new technologies or capabilities are proposed 
for high-performance aircraft, the impact, utiliza- 
tion, and costs of these technologies must be assessed 
within the context of air combat tactics and effect i ve- 
il ess. The highly complex and transient nature of air 
combat makes simulation the primary tool for per- 
forming this assessment. Both batch and real-time, 
piloted simulations can contribute to the assessment. 

Batch air combat simulations such as the ad- 
vanced air-to-air system performance model (ref. 1) 
and TAG BRAWLER (ref. 2) allow the study of air- 
craft tactics and performance in a highly controlled 
and repeatable environment. Batch air combat 
simulations consist of two fundamental elements— 
computerized maneuvering logics that generate ma- 
neuver decisions and a simulation environment in 
which maneuvering logics are developed and tested. 
Batch combat simulation programs can run large 
numbers of engagements with minimal operator in- 
tervention, which allows comprehensive sets of ini- 
tial conditions or parametric variations to be rapidly 
evaluated. Unfortunately, the minimal operator in- 
tervention inherent in batch operation slows devel- 
opment and validation of new maneuvering logics, 
which can result in relatively inflexible tactics that 
do not effectively exploit a given situation or aircraft 
capability. 

In contrast, piloted simulation provides an envi- 
ronment ideally suited for rapid tactical experimenta- 
tion and adaptation. New tactics can be investigated 


by instructing pilots to maneuver in the desired man- 
ner. Furthermore, the natural interface provided to 
the pilots encourages their participation in this devel- 
opment process and enhances their ability to assess 
the success of a given tactic. Unfortunately, because 
human pilots introduce variability, the time required 
to perform a statistically meaningful piloted air com- 
bat simulation study, combined with the availabil- 
ity and expense of the necessary facilities and pilots, 
makes a comprehensive study extremely difficult to 
perform. 

Because the strengths and weaknesses of batch 
and piloted simulations are complementary, a syner- 
gism exists when the two approaches are employed 
in concert. To fully exploit this synergy, the Langley 
Research Center is developing an integrated batch 
and piloted simulation tool known as the tactical 
guidance research and evaluation system (known as 
TiGRES in 1989 when ref. 3 was written). The Ti- 
GRES tool consists of three primary elements: an ad- 
vanced maneuvering logic that functions in real time 
and uses artificial intelligence techniques (ref. 4); a 
multidome, piloted simulation facility, the differen- 
tial maneuvering simulator (DMS, ref. a); and a 
batch simulation environment, the tactical maneu- 
vering simulator (TMS). The development and op- 
eration of the TMS and its relation to the other el- 
ements of the TiGRES tool are the focuses of this 
report. 

Unlike existing batch air combat simulation envi- 
ronments that typically use reduced-order dynamic 
models, aircraft in the TMS arc modeled using 



equations of motion, control laws, aerodynamics, and 
propulsive characteristics identical to those used in 
high-fidelity piloted simulations in the DMS. This 
commonality allows maneuvering logics developed in 
the TMS to be evaluated without modification in re- 
lation to human pilots in the DMS. The ability to test 
maneuvering logics with human pilots provides an ef- 
ficient means of validating the results of batch simu- 
lation analysis. Thus, extensive preliminary investi- 
gations of tactical maneuvering strategies, guidance 
concepts, or aircraft performance characteristics can 
be performed quickly and cheaply with the TMS. Af- 
ter t he focus of an investigation matures, a minimum 
number of piloted simulations in the DMS can con- 
firm or refine the findings of the more comprehensive 
batch analysis. 

The TMS has three basic elements. The first 
element is the model that simulates individual air- 
craft. Currently, models representative of a modern 
high-performance aircraft with and without thrust- 
vectored (TV) capability are available 1 . The second 
clement, is the tactical autopilot (TA), which enables 
maneuvering logics to command full-order dynamic 
aircraft models in bot h the TMS and DMS. The TA 
converts guidance 1 commands issued in the' form of de- 
sired angle of attack and wind-axis bank angle' into 
inputs to the inner-loop control augmentation sys- 
tem of the simulated aircraft. The 1 third (dement, 
is the TMS executive program and the synehroniza- 
tion subroutine'; these provide the capability to simu- 
late many-versus-many (MvN) air combat by running 
multiple, single-aircraft simulations in parallel. 

This report describes the capabilities and op- 
cration of the TMS. First, the background under- 
lying the development of the TMS is discussed. Next, 
the simulation environment is described. This de- 
scription details the available aircraft models, the 
TA, and the parallel implementation used to provide' 
MvN simulation. Thereafter, example engagements 
are presented to demonstrate TMS operation. The 
paper concludes with a discussion of future areas of 
research and a summary of the current work. 


Symbols and Abbreviations 


Symbols: 
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. Fy , b'z 
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lift coefficient, 

\Aft/(q x Reference area) 

mean aerodynamic chord 

force about A-, 1 -, and Z-axes. 
lb 

acceleration due to gravity, 
32.17 ft/sec 2 


h 

Jx 

Ixz 

W 


Iz 


K D a 

% 

Kla 

Kp<> 

Kp„ 

LfiF, 

L\VB 

l wk 


M 

A/ v , My ■ M/ 


M r 

in 

n 

V 


9 


9 

r 


77/ 

/ 

th 
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altitude, ft 

rolling moment of inertia, slug-ft 2 

product of inertia, slug-ft 2 

pitching moment of inertia, 
slug-ft 2 

yawing moment of inertia, 
slug-ft 2 

gain on rate of a error 

gain on rate of /i error 

gain on integral of a error 

proportional gain on a error 

proportional gain on p error 

transfer motion matrix from 
Earth to body axis 

transfer motion matrix from 
body to wind axis 

transfer motion matrix from 
Earth to wind axis 

Mach number 

moment about A7. V'-, and 
Z-axes, ft-lb 

maximum peak overshoot 

aircraft mass, slugs 

normal load factor, g units 

roll rate in body-axis system, 
dcg/scc 

pitch rate in body-axis system, 
dc^/scc 

dynamic pressure, lb/ft 2 

yaw rate in body-axis system, 
deg/sor 

Laplace operator 

body-axis components of thrust 
force, lb 

time, sec 

thrust force, lb 

velocity along A" body axis, 
ft /see 

velocity along Y body axis, ft /sec 
velocity along Z body axis, ft /sec 
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Subscripts: 

4 

a 

E 


longitudinal, lateral, and vertical 
body axes 

X-axis of inertial reference 
system 

separation along A r -axis between 
center of gravity and thrust force 
line of action 

T-axis of inertial reference 
system 

separation along T-axis between 
center of gravity and thrust force 
line of action 

Z-axis of inertial reference 
system 

separation along Z-axis between 
center of gravity and thrust force 
line of action 

angle of attack, deg 

rate of change of a, deg /sec 

angle of sideslip, deg 

flight path angle, deg 

change in azimuth angle due to 
thrust vectoring, deg 

change in elevation angle due to 
thrust vectoring, deg 

lateral stick displacement, in. 

longitudinal stick displacement, 
in. 

thrust deflection angle, deg 
body-axis pitch angle, deg 
wind-axis bank angle, deg 
rate of change of //, deg/scc 
density ratio 

maneuver plane rotation angle, 
deg 

body-axis bank angle, deg 
body-axis heading angle, deg 


aerodynamic 

aileron 

engine 


H 

horizontal stabilator 

L 

left engine 

R 

right engine 

r 

rudder 

REF 

reference 

SB 

speedbrake 

s 

stabilator 

LEF 

leading-edge flap 

TEF 

trailing-edge flap 

Abbreviations: 


ACM 

air combat maneuvering 

ACSL 

Advanced Continuous Simulation 
Language 

AML 

Adaptive Maneuvering Logic 

azimo 

engine azimuth angle as mounted 
to airframe, deg 

CAS 

control augmentation system 

DMS 

Differential Maneuvering 
Simulator 

d.o.f. 

degrees of freedom 

elcv’o 

engine elevation angle as 
mounted to airframe, deg 

MvN 

many versus many 

TA 

tactical autopilot 

TDG 

tactical decision generator 

TiGRES 

tactical guidance research and 
('valuation system 

TMS 

tactical maneuvering simulator 

TV 

thrust vectored 

lvl 

one versus one 


Background and Objectives 

During the late 1960's and 1970's, NASA funded 
the development of a computer program to provide 
an invariant or calibrated opponent for use in pi- 
loted air combat simulation studios in the newly con- 
structed DMS. (See ref, fi.) The original specifica- 
tion called for a program capable of generating tac- 
tically sound maneuver derisions and of realistically 
simulating the resulting aircraft motions for an arbi- 
trary aircraft in one-versus-one (lvl) air combat. Re- 
searchers recognized that such a program would not 
only provide an invariant opponent in the DMS, but 
could also bo used to perform rapid parametric stud- 
ies on different aircraft characteristics and to develop 
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new tactical maneuvers fur existing and proposed air- 
craft. A final requirement was for the program to 
run in real time on the computer system of the DMS 
(a Control Data 6600), which was already burdened 
with supporting the real-time, piloted simulations. 

The resulting program, the adaptive maneuvering 
logic (AML, refs. 6 and 7), distinguished itself as a 
formidable adversary against human pilots. In a real- 
time simulation with F-4 aircraft, the AML was able 
to consistently beat experienced pilots in lvl air com- 
bat maneuvering (ACM). In fact, the real-time per- 
formance of the AML is so impressive that it is used 
successfully as a training tool in several military sim- 
ulation facilities. However, to achieve real-time per- 
formance on the 19G0's vintage computer equipment 
in use in the DMS at the time, the AML has several 
key limitations that have curtailed its use except as 
an invariant opponent. These limitations have led to 
the development of TiGRES. 

Three factors severely degrade the suitability of 
the AML simulation environment for use as a re- 
search tool. First, the motion of the aircraft is 
described by a five- degree-of- freedom (d.o.f.) "p (ir_ 
formancc model, rather than a more standard ap- 
proach with six d.o.f. As described in reference 8, the 
basic idea of the performance model is to move the 
aircraft in a realistic-appearing manner during the 
transition from the current aircraft attitude to one 
that corresponds to a commanded or desired flight 
condition. In this performance model, no moments 
are calculated: therefore, no rotational differential 
equations of motion are used to model the rotational 
dynamics of the aircraft. Instead, body-axis rotation 
rates (/>, q , and r) are calculated directly as required 
to make the transition from the current body-axis 
attitude (defined by the Euler angles i\ 6 , and o) 
to the commanded attitude. The required rates are 
approximated through the following relations: 

p = (Ao - Avs'mO) /At 
q = (Atfcoso + Ate cos 0 sin o) / At > (1) 

r — (Av cos 0 cos o — A# sin e>) / At j 

where 

At’ = tVoiu — Ccur 
AO = Atom — ^cur 
A O — (Voiii — C?cur 

At = Time increment of simulation 


and the subscripts com and cur refer to command 
and current . To prevent the aircraft from rotating at 
unrealistic rates, limits are placed on the maximum 
allowable p, q< and r. If the required p, <?. or r as 
calculated from equation (1) exceeds a maximum al- 
lowable value, that value is used instead of the calcu- 
lated value. The number of d.o.f. of this performance 
model is five rather than six because the aircraft is 
always assumed to be in an attitude without sideslip, 
hence removing one d.o.f. 

The performance model greatly reduces the com- 
putation time and data storage required to simulate 
a given aircraft. The performance model also sig- 
nificantly simplifies the task of tracking commanded 
trajectories. These trajectories are characterized by 
a desired load factor n and a maneuver-plane rotation 
angle p n ,. which is defined as the angle from the neg- 
ative gravitational vertical axis — Z /.; (i.e., upward) 
to the “maneuver plane" of the aircraft. This plane 
is defined by the velocity vector of the’ aircraft and 
the net force vector (i.e.. vector sum of the gravi- 
tational. aerodynamic, and thrust forces) affecting 
the aircraft. Because by definition no unbalanced 
forces are affecting the aircraft outside the maneuver 
plane, the maneuver plane contains the trajectory of 
the aircraft. The desired n and p m can be converted 
into a corresponding body orientation for the current 
flight condition. Because the performance model al- 
lows the body rotation rates to be commanded di- 
rectly, the commanded trajectory is easily captured 
and tracked. The motion is adequate for use as an 
invariant opponent because, from the perspective of 
a pilot flying against it in a simulator, the motion 
does appear "realistic. However, to be a useful tool 
for performing analyses, the motion must be realistic 
in a physical sense rather than just appearing real- 
istic. Close-in ACM engagements consist almost en- 
tirely of transient maneuvering, and failure to model 
the dynamics of the aircraft accurately during this 
maneuvering will yield incomplete results. 

An interesting note is that the original developers 
of AML were’ well aware of the limitations of the per- 
formance model. When a sufficiently powerful com- 
puter (a Control Data Cyber l/o) became available 
in the DMS to handle a six-d.o.f. model, such a model 
was developed and compared with both the perfor- 
mance model and pilots. (See refs. 7 and 9.) The 
results of these tests showed that, although the over- 
all combat performance of the two models was simi- 
lar. significant differences (existed between the types 
of maneuvers performed bv the performance model 
and by the six-d.o.f. model. However, because the 
primary interest in AML was still on providing an 
invariant opponent, the similar combat performance 


4 



of the two models was taken as validation of the suit- 
ability of the performance model in this capacity. 
After these tests were completed, no further work 
appears to have been done with the six-d.o.f. model. 

The second deficiency of the AML simulation en- 
vironment is that it provides only for Ivl air com- 
bat simulation. Although lvl investigations are very 
useful for preliminary analysis, complications (e.g., 
cooperative tactics) of air combat that involves multi- 
ple aircraft (three or more) make multiaircraft simu- 
lations necessary to fully investigate and understand 
the effect of a given concept. The reformulation from 
an existing lvl simulation to a multiaircraft capa- 
bility would not normally be prohibitively difficult. 
However, the lack of organization and documentation 
makes this upgrade less attractive in the case of 
AML, as described in the following paragraph. 

The final deficiency of the AML simulation en- 
vironment is that the FORTRAN code that imple- 
ments the equations of motion was done in an ad hoc 
manner with various undocumented alterations and 
experiments scattered about. For example, elements 
of the maneuver decision process are implemented 
in the equations of motion routine simply because 
they were easier to implement there and may have 
increased execution speed. Having parts of the ma- 
neuver decision process scattered around in the sim- 
ulation routines not only makes following the equa- 
tions of motion more difficult, it makes tracking the 
decision process nearly impossible. Aircraft may per- 
form maneuvers in a manner that is inconsistent 
with the intended decision process because remnants 
of an earlier decision logic were “hard wired” into 
the code that implements the equations of motion. 
This convoluted code is extremely difficult to upgrade 
reliably. 

Thus, based on the need to provide a more real- 
istic air combat simulation along with the difficulty 
of upgrading the simulation environment of AML to 
meet this need, the decision was made to develop the 
TMS as a new program. The experience gained from 
working with the AML has been helpful in defining a 
set of objectives for the TMS. To support the research 
objectives of TiGRES, TMS requires the following 
features: 

1. The aircraft simulation model must be function- 
ally equivalent to models used for piloted simula- 
tion studies in the DMS. This equality will allow 
a common tactical decision generator (TDG) to 
be tested against baseline decision logics in batch 
simulations and against pilots in the DMS. Any 
differences between batch and piloted simulation 


results will be directly attributable to differences 
in maneuver strategies. 

2. Current TDG's use n and p m to characterize the 
desired trajectory. The performance model used 
by the AML allows the corresponding lift coeffi- 
cients Ci and <p to be commanded directly. Un- 
fortunately, a model that is equivalent to a piloted 
simulation model mandates the use of six-d.o.f. 
dynamics. With these higher order dynamics, the 
ability to command lift and bank angle directly 
is lost. A control system or autopilot must be 
added to the aircraft model to issue commands to 
the inner-loop control system so that the aircraft 
can capture and track the desired trajectory in 
near-minimum time. 

3. The TMS must support simulation of multiple air- 
craft. The DMS currently has hardware to simu- 
late and project three aircraft, which limits tests 
in this facility to lv2 scenarios. However, because 
future upgrades to the DMS can be anticipated, 
the structure of the TMS should accommodate 
MvN participants. 

4. The TMS must function as an independent el- 
ement, with the information flow between the 
TMS and the TDG handled in a structured and 
easily controlled fashion. This separation is in- 
tended to prevent functions of TDG's from being 
inadvertently implemented in the TMS. 

As will be shown in the following sections, the 
simulation environment described in this report 
meets these objectives. 

Tactical Maneuvering Simulator 

Functional Overview 

The TMS provides a batch simulation environ- 
ment for developing and evaluating tactical maneu- 
vering strategies. The TDG's that implement var- 
ious maneuvering strategies are tested against one 
another in varying initial conditions. The resulting 
trajectories can then be used to refine these strate- 
gies. Multiple iterations through this refinement pro- 
cess permit a globally effective maneuver strategy 
to be developed for a given aircraft. The TMS can 
also be used to evaluate the tactical implications of 
perturbations to aircraft performance or supporting 
systems. By comparing the combat performance of 
a modified aircraft (and appropriate TDG) with a 
baseline aircraft, designers can assess the effect of 
the modification. This assessment should provide an 
indication of the overall value of that modification 
in terms of an exchange ratio and the types of tacti- 
cal maneuvers and situations that favor the modified 
aircraft. 
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The TMS provides an air combat environment 
with any number of engagement participants. A 
parallel implementation structure allows individual 
aircraft simulations to be initiated or ‘"spawned” as 
needed. The number of aircraft being simulated at 
one time is limited only by the available computer 
memory and the desired computation speed of the 
simulation. Equations of motion for six d.o.f. are 
used to model the motion of each aircraft and data 
representative of a high-performance aircraft both 
with and without TV systems are available for use 
in these equations. The user is thus able to com- 
pare the performance of an enhanced agility, TV air- 
craft with that of an aircraft of conventional agility. 
The equations and data used to model the aircraft 
in the TMS are also implemented for piloted simula- 
tions in the DMS. This implementation provides the 
desired commonality between the batch and piloted 
simulation environments of TiGRES. 

The TMS has three basic elements. The first el- 
ement is the aircraft simulation model, which sim- 
ulates the motions of each participating aircraft. 
The second element is the tactical autopilot (TA), 
which controls the aircraft such that it captures and 
tracks the trajectory commanded by its correspond- 
ing TDG. The third element is the TMS executive 
program, which enables multiaircraft simulation by 
spawning individual aircraft, as needed, by over- 
seeing the engagement in a common inertial reference 
frame and by controlling communication between air- 
craft and TDG’s. These elements are described in the 
following three sections. 

Aircraft Simulation Model 

Individual aircraft are modeled with a modified 
version of an existing batch simulation model devel- 
oped at the Langley Research Center. This simula- 
tion models an F-18 aircraft with or without a hypo- 


thetical, hardware-based TV system developed by 
the Northrop Corporation. This TV system uses two 
vectoring vanes on each engine to provide thrust- 
induced pitching and yawing moments. To distin- 
guish between the aircraft equipped with the TV 
system and the basic aircraft, the basic aircraft is 
referred to as the baseline aircraft, whereas the air- 
craft with the TV system is referred to as the TV 
aircraft. The batch simulation was developed from 
the real-time simulation code for the F-18 aircraft as 
implemented in the DMS and from documentation 
obtained from the McDonnell Aircraft Company. An 
in-depth description of the batch simulation has been 
published (ref. 10), but details relevant to use in the 
TMS are presented here. 

Implementation of simulation. The com- 
puter code that implements the simulation model 
is written in the advanced continuous simulation 
language (ACSL) (ref. 11) and FORTRAN. (See 
ref. 12.) The ACSL is a simulation system with 
a special-purpose high-level language, a translator, 
and various libraries to satisfy the commands avail- 
able in the language. The ACSL simulation mod- 
els are translated into FORTRAN and linked with 
the ACSL libraries. The resulting executable pro- 
gram allows interactive user input and enables the 
generation of plots and printed outputs. The ACSL 
allows FORTRAN subroutines to be integrated into 
the simulation model. 

The simulation uses the ACSL to implement the 
dynamics of the aircraft and engines. Actuator and 
sensor models are also implemented in the ACSL. 
FORTRAN subroutines are used to calculate aero- 
dynamic forces and moments and steady-state engine 
parameters. The discrete, inner-loop, control aug- 
mentation system of the aircraft is also implemented 
primarily in FORTRAN. 


Equations of motion. The equations of motion in the ACSL simulation effectively model the flight of 
a rigid airplane over a flat, nonrotating Earth. The aircraft mass and moments of inertia are set at the start 
of a simulation and are assumed to be constant. The aircraft is considered to be symmetric about the plane 
defined by the X and Z body axes, so that the I\y and lyz products of inertia are zero and are not included 
in the equations. With these simplifications the equations take the following form: 

Translational equation 
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Kinematic relations 



1 sin 0 tan 9 cos 0 tan 9 
0 cos 0 — sin 0 

0 sin 0 sec 9 cos 0 sec 9 



( 4 ) 


Typical weights and moments of inertia used for the baseline and TV aircraft are shown in table I. Aerodynamic 
and thrust-induced forces and moments are discussed below. 


Aerodynamic forces and moments. The aerodynamic characteristics of the simulated aircraft are 
discussed in detail in references 13 and 14. Figure 1 illustrates the configuration of the aerodynamic 
surfaces and controls. Table II provides dimensional data relevant to these aerodynamic effectors. The 
aerodynamic force and moment generated by each surface or control are calculated from a large wind- 
tunnel-derived data base using table look-ups with linear interpolations. Data are stored in a non- 
dimensional form as functions of angle of attack a, angle of sideslip ft, Mach number M, the time rates 
of change of a and ft, surface deflections, and rates p, q, r. The a range is -10° to 90°, the ft range is -20° 
to 20 , and the M range is 0.20 to 2.00. Flexibility effects in the form of flex-rigid ratios and flexibility incre- 
ments are included in the data base to an altitude of 60000 ft. Actuators for all control surfaces except the 
speedbrake are modeled with a first-order lag with time constants and rate limiting, as in table II. The actuator 
responsible for moving the speedbrake is modeled as producing a constant deflection rate of 24 deg/sec. 

Engine forces and moments. Two engines rated at 16 100 lb of installed static sea level thrust are 
included in the simulated aircraft. The engine model takes inputs from the throttle and current air data 
(altitude h, dynamic pressure q, and M) to compute the force produced by the engines. For the TV aircraft, a 
and ft effects as well as thrust losses attributable to vectoring are included in the thrust computation. Given 
this information, the body-axis components of thrust for each engine are computed as 


™X R = th R cos (elev 0 + <5 e lev TVi /i) cos (azim 0 + <5azim TV , fl ) 
TH X l = th L cos (elevf) + ^ e lev TV / ) cos (azim 0 + ^ay.iiri xv / ) 
™Yr = th R cos (elev 0 + <W TV /( .) sin (azim 0 + <Wn TVJi ) 
TH Y l = ~ th L cos (elevo + <5 e ] e v TV sin (azim 0 + <5azim TV jr 
TH Zr = th R sin (elevo + <W TV /V ) sin (azim 0 + L mxv ,J 
TH Zl = th L sin (elev 0 -I- <W TV x) sin (azim 0 + <5 aziniTV / ) 


(5) 


The elevation angle of the engine is defined in the aircraft X-Z plane; positive direction is the thrust directed 
in a positive ^-direction. The azimuth angle is measured in the aircraft plane; positive direction is thrust 
directed inward toward the vehicle centerline. For the baseline aircraft, the elevation angle is 0°, azimuth angle 
is 1.98 , and the <5xv terms are 0. The TV aircraft is equipped with a TV system that has two vanes per engine 
as shown in figure 2. The change in elevation and azimuth angle produced by the TV system is defined by 

^elevxv = sin_1 (sin 48° sin 6 th ) 

^azim TV = sin_1 (cos 48° sin 6 th ) 
where S t f t is the thrust deflection angle in degrees. 

By deflecting the thrust of the two engines in a symmetric or nonsymmetric manner, a researcher can 
generate nearly pure pitching or yawing moments that are similar to those of an aerodynamic V-tail aircraft. 
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The actuators for the TV vanes are modeled as first-order transfer functions with a steady-state gain of one, 
a time constant of 1/30 sec, rate limits of 80 deg/sec, and position limits of ±30°. 

The force and moment terms in the equations of motion can now be computed as 

F Xe = T H x l + TH XM 

Fy h: = TH y L + THy M 

f Ze = th Z ' L + th Z ' R 

A/\' f ; = THz y RYe ng “ THz.L ^ eng “* Fy k 2 eilg 

My E = ~Fz E X e ng + Fx E Z e ng 

A Izy — Fy^Xeng + T H X iY en g ~~ FH X ,R^e ng , 


Control augmentation system . As a fly-by- 
wire aircraft with a full authority control augmen- 
tation system (CAS), the dynamic characteristics of 
the simulated aircraft depend heavily on the actions 
of this CAS in addition to the underlying open-loop 
dynamics described above. This CAS is documented 
in detail for the baseline aircraft in references 15 
and 16. A simulation of the “auto flap up r mode 
of the CAS defined by the version 8.3.3 produc- 
tion programmable read-only memory (PROM) set 
is in the simulation model. This auto- flap- up opera- 
tional mode of the CAS is normally engaged during 
ACM. The CAS can be divided into control loops 
about the longitudinal, lateral, and directional axes. 
The longitudinal CAS and the other two controllers 
have minimal coupling; however, the lateral and 
directional controllers are coupled through various 
interconnections and will be described together. 

The longitudinal CAS, shown in figure 3, uses the 
longitudinal stick position as the command input. 
The forward path gains are air data scheduled to 
yield a uniform initial pitch acceleration response for 
sharp stick inputs. A forward loop integrator drives 
to zero the steady-state error between the maneu- 
ver command (from longitudinal stick position) and 
the feedback variables. The CAS feedback is an air 
data scheduled blend of pitch rate, normal acceler- 
ation, and angle of attack. Pitch rate and normal 
acceleration feedbacks give improved pitch dynamic 
characteristics and load factor control in the mid- to 
high-dynamic-pressure portion of the flight envelope. 
Improved ACM flying qualities and increased stick- 
force-per-g cues in the low- to mid-dynamic-pressure 
flight regime are provided by the air-data-scheduled 
pitch rate feedback. Angle-of-attack feedback pro- 
vides additional increased stick force cues for low- 


speed, high-a ACM. Roll rate multiplied by yaw rate 
is fed to the longitudinal CAS to reduce the effects of 
inertial coupling. The longitudinal CAS also sched- 
ules the deflection of the leading- and trailing-edge 
maneuvering flaps as a function of a and air data 
to optimize performance, improve high-a character- 
istics, and provide load alleviation at elevated load 
factors. 

The lateral and directional CAS, shown in fig- 
ure 4, sums lateral stick position with roll rate feed- 
back to provide closed-loop control of the ailerons, 
differential stabilators, differential trailing-edge flaps, 
and differential leading-edge flaps. The lateral CAS 
command path consists of structural notch filters 
and air-data-scheduled gains. The gains vary with 
q, static pressure, and a to provide acceptable loop 
stability and roll response characteristics through- 
out the flight envelope. Maximum roll rate is lim- 
ited to 220 deg/sec when normal loads are less 
than 5 g and 150 deg/sec for normal loads greater 
than bg. The directional CAS uses a command sig- 
nal from the rudder pedals with stability-axis yaw 
rate (rcosa -psina) and lateral acceleration feed- 
back. The rudder pedal force transducer signal is a 
and air data scheduled to prevent a command that 
would exceed the vertical tail load limits and to elim- 
inate aircraft departures for full pedal inputs. The 
r cos a feedback component helps provide sideslip re- 
duction during moderate and high-a maneuvering 
flight. Lateral acceleration feedback aids in reduc- 
ing sideslip and provides turn coordination. Roll 
rate multiplied by pitch rate is fed to the direc- 
tional CAS to reduce the effects of inertial coupling. 
The lateral and directional controllers are coupled 
through a rolling-surface-to-rudder interconnect and 
a rudder-pedal- to- rolling-surface interconnect. The 


8 



rolling-surface-to-rudder interconnect is incorporated 
to minimize sideslip that could accompany lateral 
stick inputs. Similarly, the rudder-pedal-to-rolling- 
surface interconnect is provided to reduce sideslip 
and a excursions from rudder pedal inputs at high a. 
The interconnect is scheduled with a and is scheduled 
to zero at low a. 

The CAS used with the TV aircraft is a refined 
and extended version of the baseline CAS. This work 
was performed by the Flight Dynamics Branch at 
the Langley Research Center through extensive batch 
and piloted simulation analyses. The CAS integrates 
the TV system with the aerodynamic control sur- 
faces to significantly increase the maneuvering ca- 
pabilities of the aircraft at high a. The feedback 
structure and operation of this CAS are similar to 
those described for the baseline aircraft. The pitch 
and yaw commands from the command paths are di- 
vided, as appropriate, between the aerodynamic and 
TV controls. The pitch and yaw commands sent to 
the TV system are passed through a mixer that re- 
solves the commands into appropriate vane deflection 
commands for the TV hardware of the left and right 
engines. 

The CAS described above augments the dynam- 
ics of the bare airframe to provide stability and pre- 
dictable flying qualities that enable pilots to employ 
the aircraft in tactical engagements. For use in the 
TMS, an outer-loop control system is needed around 
the basic CAS to track trajectories commanded by 
a TDG. In a sense, this outer-loop control system 
performs the physical functions of the pilot— that 
is, it transforms the desired tactical plan or strategy 
into actual aircraft motions. This outer-loop control 
system, the TA, is described in the following section. 

Tactical Autopilot 

The TA accepts trajectory commands generated 
by a TDG and issues commands to the inner-loop 
CAS that cause the aircraft to follow the desired tra- 
jectory. Current TDG’s issue trajectory commands 
by specifying parameters that define a desired magni- 
tude and orientation for the lift force combined with 
a desired throttle and speedbrake setting. Because 
the throttle and speedbrake settings are obtained di- 
rectly, no interface is needed to capture these com- 
mands; the commands are passed directly from the 
TDG to the aircraft simulation. In contrast, the mag- 
nitude and orientation of the desired lift force cannot 
be obtained directly, which requires the development 
of the TA. 

Many different parameter pairs can be used 
to specify the desired lift vector. For a given 


flight condition, the magnitude of the lift vector can 
be specified by commands to the corresponding de- 
sired load factor to C L or to a. Similarly, the orien- 
tation can be specified by various angular references 
such as p m , 4>, or wind-axis bank angle p, which is 
defined as 

n = tan -1 ( singcoSQsin ^ + sin <p cos 8 cos 3 - cos g>cos gsin q sin P \ 
V sin# sin a + cos <2> cos 0 cos a ) 

( 8 ) 

Equation (8) is obtained from the matrices that 
transform vectors from Earth axis to body axis L BE 
(ref. 17) and body axis to wind axis L\vb (ref. 17) to 
calculate Earth axis to wind axis L\ve (ref. 18) with 
the relationship L\ve = L\vb L BE . 

For modern, high-performance aircraft, specifica- 
tion of a and p offers several advantages. First, to 
fully exploit the tactical potential of these advanced 
aircraft, the TDG must command maneuvers in 
the stall/poststall region. During those maneuvers, 
the aircraft orientation is frequently more important 
than its flight path. Because lift curve slopes are gen- 
erally shallow and variable in the stall/poststall re- 
gion, orientation relative to velocity vector is poorly 
defined by load factor and C E . In contrast, a remains 
an effective command variable in the stall/poststall 
region. Second, an awareness of a is ensured in the 
TDG. Because the current and future maneuvering 
potential of an aircraft is largely a function of a, 
this awareness is imperative to the formulation of 
effective maneuver decisions and strategies. Third, 
p directly specifies the desired orientation of the lift 
vector, thereby eliminating the need to calculate the 
corresponding body-axis bank angle while ensuring 
that the vector is oriented as intended. 

The TA thus is an all-attitude, outer-loop con- 
trol system to capture and track o and p as com- 
manded by a TDG. Coordinated flight (defined as 
flight with 0 = 0) is assumed desirable at all times. 

A block diagram of the complete TDG TA aircraft 
system is shown in figures 5 and 6. The TA described 
in this paper represents an initial design and allows 
current TDG’s, intended to operate with five-d.o.f. 
performance models, to interface with and effectively 
command full six-d.o.f. models. The TA enables this 
interface with minimal modifications to these exist- 
ing TDG’s. As experience is gained from these initial 
efforts, the design of the TA can be refined as per- 
formance requirements and even desired command 
variables become better defined. For instance, full 
exploitation of the nose-pointing capability of the 
simulated aircraft may make 3 = 0 not desired at all 
times. 


9 



The task performed by the TA is similar to the 
function of the control system developed for the six- 
d.o.f. model test in the AML. This control system, 
which is described in reference 9, allowed the guid- 
ance logic of the AML to effectively command a six- 
d.o.f. simulation of an F-4 aircraft. Because of this 
success and the similarity to the current application, 
reference 9 has been a guide during the development 
of the TA. The design and development of the TA is 
described in detail in reference 19 and is summarized 
herein. Although the TA is described in this report in 
the context of the TMS, its use is also required in the 
DMS. The incorporation of the TA into the piloted 
simulation model of the DMS permits the TDG s to 
command this simulation in an identical manner to 
the batch simulation. 

The TA is divided into two channels - a longitu- 
dinal command system that uses longitudinal stick 
inputs to capture and track commanded a and a lat- 
eral command system that uses lateral stick inputs to 
capture and track the commanded p. A directional 
controller is not included in the TA because the inner- 
loop CAS already attempts to maintain zero sideslip, 
unless commanded otherwise by the rudder pedal in- 
puts. Piloted simulations have shown that the wind- 
axis rolling performance of the baseline aircraft can 
be improved slightly at a > 25° by rudder pedal in- 
puts. (See ref. 20.) This performance is not being 
exploited by the current TA. 

The longitudinal command system uses a 
proportional-integral-derivative (PID) structure with 
a feedback, as shown in figure 6(a). The lateral 
command system uses a proportional-derivative (PD) 
structure with p feedback, as shown in figure 6(b). 
The values of a, the rate of change of a (d), p, and 
the rate of change of p (p) are assumed to be avail- 
able without error, so no additional compensation to 
account for sensor noise or dynamics is included in 
the TA. Also, no attempt is made to model the cog- 
nitive and neuromuscular delays or limitations that 
arc inherent in a human pilot. Thus, as implemented, 
the TA represents an idealized controller. 

The gains for the command systems were de- 
termined through a combination of linear analysis 
and evaluation of the full nonlinear system response 
to step commands and representative command se- 
quences. To obtain good performance throughout 
the ACM envelope of the simulated aircraft, the three 
gains of the longitudinal command system {Kp a , 
Kp n , and K In ) are scheduled as a function of q. 
In addition, Kp a is also scheduled <is a function of 
density ratio p/p G to compensate for changes in aero- 
dynamic damping with altitude. Good performance 
across the ACM envelope is achieved by the lateral 


command system by the scheduling of its two gains 
Kp M and Kp^ with a. 

To achieve time-optimal control of a system with 
limited control authority, generally the maximum 
available control authority must be used at all times. 
(See ref. 21.) Because the TA should capture com- 
mands in minimal or near-minimal time, the gains of 
the command systems have been selected such that 
the commanded stick positions are frequently near 
saturation for small command changes and saturated 
for moderate and large changes. This saturation does 
not cause significant difficulties for the lateral com- 
mand system. Gains Kp ^ and Kp^ are selected such 
that the lateral stick input becomes unsaturated with 
sufficient control authority remaining for the linear 
controller to capture the desired p with acceptable 
levels of overshoot. Saturation can cause problems 
with the longitudinal control system unless the ac- 
tion of the integral clement is restricted to prevent 
integrator windup. If the gain on the integral cle- 
ment is adjusted such that good response is achieved 
for small command changes, large overshoots are ob- 
tained for moderate and large changes. During these 
changes, the maximum rate is quickly reached at 
which a can be increased (or decreased). Because of 
this nonlinear, rate-limited performance, the longitu- 
dinal stick command from the integral control action 
can reach very high levels during the initial response. 
The integral of the a error decreases only after the 
desired a is exceeded, so large overshoots can result. 
To prevent this windup, the calculation of the inte- 
gral of the a. error is suspended when the sum of the 
longitudinal stick commands from the proportional 
and rate elements causes saturation. This suspen- 
sion is bypassed if the current integral command is 
in opposition to the direction of saturation. This by- 
pass is necessary to efficiently respond to command 
changes that involve a sign change in a error. 

During evaluations of system response to coupled, 
large-amplitude a and p commands, the baseline air- 
craft was discovered to be prone to departures at rel- 
atively low a when full or nearly full lateral stick 
inputs are used and when the longitudinal stick in- 
put is aggressively increased to maintain constant a. 
As shown in figure 7, the departure results because (5 
builds to excessively high levels as the rudders satu- 
rate against their deflection limits. This departure 
results when the inner-loop CAS allows the air- 
craft to obtain a roll rate beyond its ability to re- 
main coordinated. As the departure represents a 
potentially dangerous flight characteristic, the phe- 
nomenon was investigated further in piloted simu- 
lation with the DMS. A similar, but less violent 
response was reproduced in the piloted simulation. 
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The (3 departure occurred only after the aircraft had 
rolled through 360°. In tactical maneuvering, full lat- 
eral stick will not likely be maintained much beyond 
a 180° roll; thus, this performance is unlikely dur- 
ing normal operations. The difference in departure 
characteristics observed in the TMS and DMS may 
be caused by the abrupt control commands issued 
by the TA versus those of a human pilot. To prevent 
the baseline aircraft from departing while under the 
control of the TA, the allowable stick input must be 
limited in the affected a range. For a < 15°, the in- 
put is limited to 85 percent of the maximum lateral 
stick travel. For a > 15°, the limit is relaxed in a 
linear fashion until full travel is available at a = 20°. 

One difficulty in developing a system such as the 
TA is the determination of suitable criteria with 
which to measure the acceptability of the final de- 
sign. Traditional performance specifications such as 
frequency and damping are inappropriate because of 
the large-amplitude, coupled maneuvers performed 
by the TA. Criteria that reflect the nonlinearities of 
the task must be used to assess TA performance. 
The intent of these criteria is to ensure that the 
TA can capture and track commands from the TDG 
without adversely biasing the tactical performance 
of the TOO TA aircraft system. This tactical per- 
formance is dependent on the combined interactions 
of all three components, so the response of the TA 
aircraft system should be c haracterized in relation to 
some functional benchmark. Because the only previ- 
ous controllers to demonstrate mastery of the simu- 
lated aircraft in ACM are human pilots, the perfor- 
mance of pilots with representative maneuvers can 
provide a benchmark for TA performance. 

Tables III and IV show the minimum and average 
time required for a series of experienced pilots to 
perform largo-amplitude, decoupled o and p captures 
in the baseline’ and TV aircraft, as simulated in the 
DMS. Also shown in the tables is the time required 
by the TA to perform the same captures. Time 
histories for these TA maneuvers are presented in 
figures 8 and 9. All runs start from \<j level Might 
and end when the desired cv or p is captured within 
the specified tolerance. The' tables show that for 
all but two of the tasks, the TA required less time 
than did t. lie pilots. The TA is probably able to 
consistently perform the desired maneuvers in less 
time than the human pilots because it can respond 
instantly to the current situation. In the two tasks 
in which the TA did not outperform the pilots, the 
performance differences arc’ small. 

For the 90° roll maneuver at cv = 10° with the 
TV aircraft, the TA takes 0.06 sec longer than the 
minimum piloted time. This increase is probably 


tactically insignificant and may be attributable to 
a variations during the maneuver. Data recorded 
during the maneuver show that the pilot allowed 
the a to fall to 7.2° during the maneuver: the TA 
minimum a was 8.5°. 

For the capture task at M = 0.60 and a = 40° 
with the baseline aircraft, the TA was unable to pre- 
vent the initial overshoot from exceeding the desired 
±2.0° capture tolerance. This overshoot increased 
the capture time of the TA for the original capture 
tolerance beyond the minimum piloted time. The 
initial TA overshoot was 0.44° beyond the desired 
capture tolerance. As this overshoot only slightly 
exceeds the desired capture tolerance, the tactical 
performance should not be significantly affected. Be- 
cause attempts to improve the response at this one 
condition resulted in an overall decrease in system 
performance, the decision was made to accept the 
nominal response of the system. The time listed in 
table III represents the performance of the TA with 
the capture criteria relaxed to 2.44°. 

Also shown in the table is the maximum peak 
overshoot hip for the TA captures. Burgin and 
Eggleston (ref. 9) recommend that for good tactical 
performance, M P for decoupled inputs should be 
limited to 5° in pitch and 20° in roll, regardless of 
the amplitude of the input. For all the captures, the 
TA is below these recommended limits. 

The capture tasks shown in tables III and IV 
measure performance for single-axis, step inputs. In 
ACM, the TA will be expected to respond to se- 
quences of simultaneous o and p commands. The 
responses of the TA to a representative command 
sequence are shown in figures 10 and 1 1 for the base- 
line and TV aircraft, respectively. These command 
sequences were obtained by discretizing, at 1-sec in- 
tervals, continuous a and p time histories recorded 
during piloted ACM engagements. This discretiza- 
tion was performed to obtain command sequences 
that are representative of the command update rate 
of a typical I DC. Because these command sequences 
were obtained from actual trajectories, the sequences 
should be reasonably close to the capabilities of 
the TA-cont rolled aircraft and representative of a 
tactically realistic command sequence. 

The TA appears to follow both sequences with 
sufficient accuracy to effectively implement realistic 
maneuver sequences. As shown in figures 10 and 1 L 
the ability of the TA to capture and maintain o and // 
is only slightly reduced bv the coupled command se- 
quences. However, an absolute, operational assess- 
ment of TA effectiveness cannot be performed until 
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the system is interfaced with an appropriate TDG 
and tested against human pilots in the DMS. 

Multiple Aircraft Simulation and TMS 

Executive Program 

The TMS uses a novel parallel implementation 
technique to provide multiaircraft simulations. Most 
batch multiaircraft simulation environments are im- 
plemented as a single large process. A main program 
calls various subroutines to implement the engage- 
ment participants. The researcher can create addi- 
tional participants by duplicating the requisite sub- 
routines, by renaming variables and common blocks 
as necessary to avoid memory conflicts, and by up- 
dating the calling sequence of the program. The TMS 
exploits parallel processing libraries provided by the 
Digital Equipment Corporation VAX/VMS 5.0® op- 
erating system (ref. 22) to implement simulation par- 
ticipants as independent processes that communi- 
cate with and are synchronized by a master process 
through a shared block of memory. This implemen- 
tation allows a single copy of the simulation program 
to run concurrently as needed to simulate the indi- 
vidual engagement participants. Because they are 
run as independent processes, memory conflicts are 
avoided without the need to manually modify each 
participant. The number of concurrent copies of the 
simulation that can be executed simultaneously is 
limited only by available computer memory and the 
desired execution speed of the simulation. Of course, 
an appropriate TDG would be needed to command 
the aircraft. 

In addition to simplifying the simulation of mul- 
tiple aircraft, this parallel implementation offers sev- 
eral other key advantages compared with conven- 
tional methods. Because all aircraft are simulated 
by the same program, corrections or updates to this 
model need only be performed once, which eases con- 
figuration control issues. With a conventional imple- 
mentation, these changes must be repeated in each 
duplicated subroutine. This need to repeat changes 
is frequently a source of difficulty, as the odds of a 
programming error increase with each repetition. As 
will be shown, the current parallel implementation al- 
lows different simulation models to be incorporated 
into the TMS and be intermixed with the current 
aircraft simulation model with only the addition of 
a standard subroutine. Thus, simulations of differ- 
ent aircraft types can easily be added to the TMS, 
which allows comparisons of the tactical performance 
of different types of aircraft. Simulations that may 
be added to the TMS are not restricted to aircraft; 
for example, high-fidelity missile simulations could 
also be implemented in a similar fashion. Finally, al- 


though not investigated in this study, parallel imple- 
mentation should allow individual simulations to be 
distributed on multiple, networked computers. Thus, 
if the number of simulation participants grows be- 
yond the capacity of a single computer, the ability to 
use distributed processing on an existing computer 
network may obviate the need to purchase a more 
powerful computer. 

The concurrent parallel implementation provides 
the above-mentioned benefits, but a control mecha- 
nism is needed to synchronize the otherwise indepen- 
dently executing simulations. This synchronization 
is required so that the simulations remain together 
on the same time step. Because the simulations ex- 
ecute as independent processes on a given computer 
(or computers), the order and length of time in which 
the computer operates on each process are functions 
of other jobs on the machine and are essentially inde- 
terminate. Thus, without some type of control mech- 
anism, the simulations would progress at different 
rates. 

The TMS uses barrier synchronization to control 
the progress of individual simulations. Barrier syn- 
chronization involves the use of barrier statements 
that suspend execution of individual processes at 
a specified point until all relevant processes have 
reached their respective barriers. After all processes 
have reached the barrier statements, the processes 
are allowed to continue execution. Barriers are used 
in the TMS to suspend the execution of the aircraft 
simulations at the end of the current time step or 
simulation frame. The simulations are allowed to 
proceed only after all simulations have reached the 
end of the current time step. 

The key elements of the parallel implementation 
used by the TMS are a FORTRAN executive pro- 
gram and a FORTRAN subroutine that was added 
to the aircraft simulation model to communicate 
with the executive program and to enable the exec- 
utive to synchronize the concurrently executing sim- 
ulations. The executive program is a master pro- 
cess that initializes the individual simulation models 
and supervises their operation. The executive pro- 
gram also handles communication with the TDG’s 
and passes information to and from the TDG’s by 
means of subroutine calls. Because all communica- 
tion between a TDG and its corresponding aircraft 
must pass through the executive program, the flow 
of information can be closely monitored and con- 
trolled. The final function of the executive program 
is to track and “score” the engagement in a com- 
mon reference frame. The executive program uses 
data returned from the simulations to determine the 
current relative geometry between aircraft. These 
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relative geometries are used to score the engage- 
ment by calculating the probability that each air- 
craft will successfully fire a weapon at opposing air- 
craft. This probability of kill P k is currently based on 
very simple models of the firing envelopes of a mod- 
ern, all-aspect, air-to-air missile and a high-velocity 
gun. The operational interaction between the ex- 
ecutive program, the aircraft simulation model, and 
the TDG is shown graphically in figure 12 and is 
described below. 

The TMS executive program consists of two pri- 
mary sections of code. The first section is presented 
in simplified form in appendix A and sets up the 
area of shared memory used to communicate with 
the other processes. This memory is contained in 
the common block SHARED .DATA. This common 
block is analogous to a standard FORTRAN com- 
mon block, but rather than being shared among 
subroutines of a single process, this common block 
can be shared by concurrently executing processes. 
Next, a do-loop is used to initialize each simula- 
tion participant. The command files executed by 
the LIBSSPAWN command assign unique input and 
output files to each aircraft simulation. Each time 
the PPLSSPAWN command is performed, the exe- 
cutable code of the simulation model (F18XX.EXE) 
is initialized as a new process. The command 
PPL$WAIT_AT.BARRIER( BARRIER JNT) keeps 
the simulation from proceeding prematurely and 
causing difficulties during the assignment of input 
and output files. A corresponding barrier is in the ini- 
tialization code of the simulation model. At the com- 
pletion of this first section of code, the simulations 
have been initialized and are waiting to continue 
execution at time zero. 

The second section of the executive program 
maintains the synchronization of the simulations, 
scores the engagements, and calls the TDG’s at 
each time step. This second section of code inter- 
acts with the previously mentioned communication 
and synchronization subroutine. This subroutine, 
shown in appendix B, is implemented in the sim- 
ulation model as the last routine to be executed. 
Just before the individual simulations reach the bar- 
rier BARRIER-DATA, the data shared with the ex- 
ecutive program are updated to the current time 
step. These data include the current attitude, po- 
sition, velocity, rotation rates, control positions, and 
thrust of the aircraft. The data from a specific air- 
craft can be identified by MYJNDEX. As each pro- 
cess is spawned, the operating system assigns it a 
unique integer index that can be retrieved by the 
command PPLSGET .INDEX. After all the simula- 
tions have reached BARRIER-DATA, the executive 


program is allowed to proceed to the relative geom- 
etry and P k calculations and to communicate with 
the TDG’s. The TDG’s return updated maneuver 
commands in the form of desired a, p, throttle posi- 
tions, and speedbrake settings. During this interval, 
the simulations are held at BARRIER.CMD. When 
the executive program completes this communication 
and reaches BARRIER.CMD, the simulations are al- 
lowed to proceed and receive the updated maneu- 
ver commands through the shared common block. 
It is important to recognize that the communica- 
tion and synchronization subroutine could be incor- 
porated into most ACSL or FORTRAN simulations, 
so that many different simulations can be added and 
mixed in the TMS with minimal effort . Of course, be- 
cause the TA is aircraft dependent, it would require 
retuning or redesigning to support other aircraft. 

The following section demonstrates the capabili- 
ties of the TMS through two sample engagements. 

Demonstration of Tactical Maneuvering 
Simulator 

The operation of the TMS is demonstrated by 
two example engagements. The first example demon- 
strates TMS simulation and synchronization of four 
aircraft. The second example demonstrates a lvl 
engagement between a drone aircraft that follows a 
predefined command sequence and an actively guided 
aircraft. 

Simulation of Four Aircraft 

The parallel implementation in the TMS provides 
an efficient and flexible environment for simulating 
multiple aircraft. However, because a parallel im- 
plementation introduces the possibility of synchro- 
nization problems not found in serial programming, 
the barrier structure must be specifically tested to 
ensure that no unanticipated conflicts or problems 
occur. The following example is designed to demon- 
strate the simulation of four aircraft and to check for 
proper synchronization. 

A simple maneuvering logic was developed to 
cause an aircraft that flies down the X E -axis in a 
negative direction to perform a vertical reversal ma- 
neuver, shown in figure 13. This reversal consists 
of a half-loop followed by a 180° roll to return to 
upright, level flight. The maneuvering logic divides 
the reversal into four phases. In the first phase, the 
aircraft maintains 1 g trimmed flight. In this ex- 
ample, the trim conditions are M = 0.90 at an al- 
titude of 10000 ft. The second phase of the maneu- 
ver begins when the aircraft passes over the Y^-axis. 
During this second phase, a is commanded to 10° 
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while fi = 0°. The throttle is also maintained at its 
trimmed position during this initial pull-up. The 
third phase of the maneuver begins when the flight 
path angle 7 passes through 90°. At this point in 
the trajectory, the actual fi flips from 0° to 180°. To 
maintain the aircraft in the desired pull-up, the com- 
manded is also flipped to 180°. During this phase, 
the aircraft is flying in an inverted orientation rela- 
tive to the inertial reference system. To circularize 
the trajectory, the commanded a is reduced to 6 ° and 
the throttle is increased to full afterburner. The third 
phase of the maneuver begins when 7 passes back 
through 15°. At this point, both o and ft are com- 
manded to 0°. These commands cause the aircraft 
to roll 180° from an inverted to an upright orienta- 
tion relative to the inertial system. The final phase 
of the maneuver begins when this rolling command 
is completed. To resume approximately level flight, 
q is commanded to 3 ° and the throttle is reduced to 
just above its original trimmed position. The bank 
angle is commanded as needed to remove any lateral 
offset during the 180° roll. 

The input or trim file used to provide the initial 
conditions for the simulated aircraft at the start of 
this maneuver is shown in appendix C. This trim 
file is read by the simulation and specifies the initial 
aircraft characteristics and flight condition. The file 
allows the user to vary inqrtial properties, select vari- 
ous modeling options, and specify the initial position 
and flight conditions. As shown in appendix C. the 
aircraft in this example is initialized with the inertial 
properties of the baseline aircraft and the modeling 
options are set to duplicate the DMS real-time sim- 
ulation. The flight condition is specified as straight 
and level flight at M = 0.90 and h = 10 000 ft. The 
initial position for the aircraft is set to A/,; = 5000 ft, 
Yg — 0, and c- = 180°. 

The ability of the TMS to simulate and synchro- 
nize multiple aircraft is demonstrated by using the 
maneuver commands for this one aircraft to com- 
mand three additional aircraft, starting in symmetry 
on the A g- and Y/r-axes and converging toward the 
Xg.Yg origin. As the original aircraft performs the 
reversal, the 0 . /*, and throttle commands are echoed 
to the new aircraft.. The original aircraft is in a po- 
sition that would be analogous to the flight loader 
of an aerobatic demonstration team calling out com- 
mands for the other team members to follow with- 
out question. The initial conditions and execution of 
this maneuver are such that if proper synchroniza- 
tion is maintained, the aircraft will simultaneously 
pass over the Xg,\g origin at the top and bottom 
of the reversal maneuver. 


The TMS was configured to spawn four copies of 
the aircraft simulation. Trim files identical to the 
one shown in appendix C with the exception of the 
initial Xg> Y'g, and ip were created for the three ad- 
ditional aircraft. The values of Xg y Yg. and ip of 
these trim files were set to provide the desired start- 
ing symmetry about the Xg , Yg origin. As the origi- 
nal aircraft performed its reversal, its commanded a, 

/i, and throttle positions were passed through the 
TMS executive to the other three simulations. Thus, 
if synchronization is maintained in the TMS, the re- 
sulting trajectories should remain symmetrical about 
the origin and because of the geometry of the maneu- 
ver, the four aircraft should “’collide 1 ’ at the top and 
bottom of the maneuver. Figure 14 shows the trajec- 
tories of the aircraft during the maneuver from vari- 
ous perspectives. As can be seen from that figure, the 
reversals are completed with complete symmetry and 
expected intersections, and demonstrate that correct 
synchronization is maintained. 

One-Versus-One Engagement 

The second example engagement demonstrates 
a lvl dogfight between a drone aircraft in a pre- 
defined, open-loop command sequence and an air- 
craft actively guided by a simple TDG. The objec- 
tive of this example is to demonstrate the operation 
of the TMS with a fully active TDG. 

The TDG commands a and (.1 to cause the flight 
path of the guided aircraft to intersect a. predicted 
future position of the drone aircraft. This predicted 
future position is obtained by extrapolation along a 
second-order curve fit to the past three recorded po- 
sitions of the drone aircraft. The TDG then deter- 
mines the' maneuver plant' and load factor required to 
intercept that position given the current state of the 
guided aircraft. The required maneuver plane and 
load factor are converted into a required a and //. If 
the required load factor is outside the aerodynamic 
or structural capabilities of the aircraft, the o that 
corresponds to maximum available or allowable lift 
is commanded. In addition, if the commanded // dif- 
fers from the current g by more than 45° and the 
commanded o is greater than 15°. the o command 
is reduced to 15° to expedite the execution of the 
rolling maneuver. This reduction was heuristieally 
selected and does not necessarily reflect an optimal 
maneuvering strategy. 

The engagement, between the two aircraft is 
shown in figure 15 from various perspectives. The 
engagement starts with both aircraft trimmed in 1 g 
level flight at h = 10 000 ft and M = 0.90. Both air- 
craft start from opposite headings with a longitu- 
dinal separation of 10 000 ft and a lateral offset of 
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1000 ft. The drone aircraft is initially commanded 
to maintain p = 0° and to increase a slightly over 
the trim value. The throttle of the drone aircraft is 
advanced into the afterburner region. These com- 
mands are maintained during the first 10 sec of the 
engagement. After the initial merge, the guided air- 
craft responds by performing an oblique, pitch-back 
maneuver to reverse its heading back toward the 
drone aircraft. After this initial period, the drone 
is commanded to increase a to 28° and to alternate 
p between ±90°, switching every 10 sec. The re- 
sulting motion is a descending spiral trajectory. In 
response to these maneuvers, the guided aircraft re- 
verses its heading again and effectively tracks the 
drone down the descending spiral. Time histories 
of commanded n versus actual o and commanded p 
versus actual p for the guided aircraft are shown in 
figure 16. These time histories demonstrate that the 
TA-controlled aircraft can closely track the TDG- 
generated guidance commands. 

These two examples have demonstrated the op- 
eration of the completed TMS. The following sec- 
tions describe potential future research activities and 
summarize the accomplishments of the current work. 

Future Research Activities 

Future research options include development of 
additional aircraft simulation models, incorporation 
of human physiological factors into the design of the 
TA, and the addition of an interactive user inter- 
face to allow the TMS to function as a tactical 
workstation. 

Because the parallel implementation technique al- 
lows aircraft simulations to be added to the TMS 
with minimal effort, numerous existing simulations 
could be added to the environment, thereby provid- 
ing the user with a catalog of aircraft types. An 
interesting model to include in this selection would 
be an unmanned aircraft, designed without the phys- 
iological and safety constraints imposed by a human 
pilot. A very illuminating test could be run that com- 
pares the performance of this type of aircraft, flown 
by a TDG, with conventional piloted aircraft. Use of 
this unmanned aircraft as an “automated wingman” 
to support conventional piloted aircraft could also be 
investigated. 

The basis for the current TA was the assumption 
that the inner-loop control system of the aircraft pro- 
vides desirable handling qualities. This assumption 
could be tested further by incorporating elements of 
pilot modeling into the TA. The field of pilot model- 
ing is an attempt to quantify the controlling actions 
of a pilot through appropriate transfer functions. 


Terms are incorporated into these transfer functions 
to reflect the physical capabilities and limitations of 
a typical pilot. Existing theory is limited largely 
to control of a single axis for small-amplitude track- 
ing tasks and significant research would be required 
to extend this theory throughout the TA operating 
range. If successful, the TMS could provide an initial 
assessment of the combat effectiveness of preliminary 
or proposed aircraft designs as flown by a typical pilot 
in tactical engagements. This initial assessment has 
several advantages: it could be performed quickly, 
it would be inexpensive, it would reduce the need 
for piloted simulation, and it would allow designers 
to make more informed decisions during the design 
process. 

The TMS currently depends on TDG’s to gener- 
ate trajectory commands for the simulated aircraft. 
However, the TMS could be easily modified by the 
addition of an interactive user interface to receive 
commands from human operators for some or all air- 
craft. The TMS could thus be used as a tactical 
workstation, allowing pilots and tacticians to explore 
maneuvering strategies in low-cost, nonreal-time sim- 
ulations. The ability to bring the human clement 
into ACM studies during the batch simulation phase 
should significantly reduce the time required to val- 
idate results in real-time, piloted simulations. To 
maintain the situational awareness necessary to de- 
velop effective maneuver strategies, these operators 
will need a large quantity of data, which can probably 
be conveyed most efficiently by a graphical interface. 
Ideally, this interface would allow pilots who are un- 
familiar with the system to intuitively and effectively 
command simulated aircraft after a brief instruction 
period. 

Concluding Remarks 

The development and operation of a batch air 
combat simulation environment known as the tactical 
maneuvering simulator (TMS) have been presented. 
The TMS is a tool for developing and evaluating 
tactical maneuvering logics. The environment can 
also be used to evaluate the tactical implications of 
perturbations to aircraft performance and supporting 
systems. 

The TMS was developed from an existing batch 
simulation of a modern, high-performance aircraft, 
with and without thrust vectoring. This batch sim- 
ulation uses six-degree-of-freedoin (d.o.f.) equations 
of motion, aerodynamics, propulsive characteristics, 
and control laws equivalent to those in high-fidelity 
piloted simulation. 

An outer-loop control system, the tactical auto- 
pilot (TA), was developed to allow existing guidance 
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logics intended for use with a reduced-order aircraft 
model to command the six-d.o.f. aircraft model with 
minimal modification. The TA uses longitudinal and 
lateral stick inputs to capture angle of attack and 
wind-axis bank angle as commanded by the guid- 
ance logic. The performance of the TA was demon- 
strated by comparison of the time required for it 
to capture decoupled angle-of-attack and bank-angle 
commands with the time required by human pilots 
for the same commands. The TA performed as well 
as or better than the pilots for nearly all the com- 
mands investigated. The ability of the TA to track 
realistic command sequences of angle of attack and 
bank angle was demonstrated on sequences gener- 
ated from piloted air combat simulations. The TA 
was shown to effectively track these representative 
command sequences. 

To provide for the simulation of air combat 
with multiple participants, a parallel implementa- 
tion scheme was developed from the parallel pro- 
cessing libraries provided by the Digital Equipment 
Corporation VAX/ VMS 5.0® operating system. This 
parallel implementation allows the TMS to simulate 
air combat with any number of engagement partici- 
pants; in fact, the maximum number is limited only 
by the available computer resources. The parallel 
implementation also simplifies software maintenance 
and allows new simulations to be easily added to the 
environment. 

The capabilities of the TMS were demonstrated 
with two example engagements. The first engage- 
ment demonstrated TMS ability to simulate four 
aircraft; the second demonstrated TMS ability to 
interact with an active guidance logic. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
April 27, 1993 
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Appendix A 

TMS Executive Program 

The TMS executive program is presented here in simplified form. It is shown dimensioned for up to four 
aircraft. The function of this routine is to initialize engagement participants and oversee the engagement in a 
common reference frame. 

PROGRAM TMS _EXEC 
C 

C EXTERNAL DEFINITIONS 
C 

INTEGERS PPL$SPAWN , LIB$SPAWN, PPL$INITIALIZE 
INTEGER* 4 PPL$CREATE -BARRIER, PPL$WAIT-AT -BARRIER 
INTEGER*4 PPL$CREATE .SHARED -MEMORY , LIB$PUT .OUTPUT 
C 

C LOCAL DATA 
C 

INTEGER* 4 LENADRC2) , STATUS 
INTEGER*4 ONE .PAGE 
PARAMETER (0NE_PAGE = 512) 

C 

REAL RANGE (4, 4) , RANGE-RATE (4, 4) 

REAL LOS (4, 4) , AZIMUTH(4,4) , DEVIATION (4, 4) , ANGLE_0FF(4,4) 

REAL MIS_PK(4,4) , GUN_PK(4,4) 

C 

C DATA FOR SHARING 
C 

BYTE FRONT -GUARD (ONE-PAGE) 

INTEGER COPIES 

REAL AIRSPEED (4) , ALPHA (4) , BANKWND(4) , BETA(4) 

REAL DIRC0S(4,9) , EULER(4,3), GAMMA (4) 

REAL GL0AD(4) , MCH(4) , P0SITI0N(4,6) 

REAL qUAT(4 ,4) , R0TRATES(4,6) , SPDBRAKE(4) 

REAL STKRUD(4,3) , TIME(4) , THRUST(4) 

C 

REAL C0M_ALPHA(4) , C0M_BANK(4) , C0M_SPDBRK(4) , C0M_THRUST(4) 

BYTE REAR_GUARD (0NE_PAGE) 

C 

C PUT SHARED DATA IN TO COMMON BLOCK 
C 

COMMON/SHARED _DATA/FRONT_GUARD , 

1 COPIES, 

1 AIRSPEED, ALPHA, BANKWND , BETA, 

1 DIRCOS, EULER, GAMMA, 

1 GLOAD, MCH, POSITION, 

1 QUAT, ROTRATES, SPDBRAKE, 

1 STKRUD , TIME, THRUST, 

C 

1 C0M.ALPHA , C0M.BANK , 

1 COM-SPDBRK, COM_THRUST, 

1 REAR_GUARD 
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c 

CHARACTER* 8 PLANE (4) 

DATA PLANE/ ’ QPLANE1 ’ , * QPLANE2 ’ , ’ 0PLANE3 ’ , ’ 0PLANE4 ’ / 

C 

C MAP SHARED ADDRESS SPACE 
C 

LENADR(l) = V.LOC (REAR-GUARD) + ONE .PAGE - V.LOC (FRONT .GUARD) 
LENADR(2) = V.LOC (FRONT .GUARD) 

PRINT * , ’ PEND LENADR’ ,LENADR(1) ,LENADR(2) 

STATUS = PPL$CREATE .SHARED .MEMORY (SHARED .DATA , LENADR) 

PRINT * , ’PEND LENADR’ ,LENADR(1) ,LENADR(2) 

C 

C LOOP TO CREATE AIRCRAFT 
C 

STATUS = PPL$CREATE_BARRIER(BARRIER.INT, ’ BARRIER. I NT’ , V.REF(2)) 
PRINT *, ’INPUT NUMBER OF AIRCRAFT (1-4). ’ 

READ (5,11) COPIES 
11 FORMAT (I 2) 

DO 99 I = 1, COPIES 

IF (I.EQ.l) STATUS = LIB$SPAWN ( ’ 0PLANE1 ’ ) 

IF (I.EQ.2) STATUS = LIB$SPAWN ( ’ 0PLANE2 ’ ) 

IF (I. EC). 3) STATUS = LIBSSPAWN ( ’ QPLANE3 ’ ) 

IF (I.EQ.4) STATUS = LIB$SPAWN ( ’ 0PLANE4 ’ ) 

N=1 

STATUS = PPL$SPAWN(N, ’ [KHG.SIM.XTMS.F18XXDF18XX.EXE’) 

STATUS = PPLSWAIT .AT .BARRIER (BARRIER.INT) 

99 CONTINUE 
C 

STATUS = PPLSCREATE _BARRIER(BARRIER_DATA, ’ BARRIER J)ATA’ , 

'/.REF (COPIES+1) ) 

STATUS = PPL$CREATE_BARRIER(BARRIER_CMD, ’BARRIER.CMD’ , 

‘/.REF (COPIES+1)) 

TSTP =90.0 
I STEP = TSTP * 32 
INITIAL = 1 
C 

C OPERATE LOOP 
C 

DO 101 I = 0 , I STEP 
C 

STATUS = PPL$WAIT_AT .BARRIER (BARRIER .DATA) 

C 

CALL PKILL( RANGE, 

RANGE .RATE, 

LOS, 

AZIMUTH, 

DEVIATION, 

ANGLE.OFF, 

MISJ’K, 

GUN _PK ) 
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CALL TMS.l (desired input, A_C0M1, B_C0M1 , THR0T-C0M1, SPDBRK.COM1) 
COM-ALPHA(l) = A.COM1 
COM_BANK(l) = B.COM1 
COM_THRUST(l) = THROT_COMl 
COM-SPDBRK ( 1 ) = SPDBRK-COM1 
C 

CALL TMS.2 (desired input, A .COM2, B.C0M2, THR0T-C0M2, SPDBRK_C0M2) 
C0M_ALPHA(2) = A.C0M2 
C0M_BANK(2) = B.C0M2 
COM_THRUST (2) = THR0T_C0M2 
COM_SPDBRK (2) = SPDBRK-C0M2 
C 

CALL TMS_3 (desired input, A.C0M3, B.C0M3, THR0T-C0M3 , SPDBRK.C0M3) 
C0M_ALPHA(3) = A_C0M3 
C0M_BANK(3) = B.C0M3 
COM-THRUST (3) = THR0T.C0M3 
COM_SPDBRK (3) = SPDBRK.C0M3 
C 

CALL TMS_4(desired input, A.COM4, B.C0M4, THROT_COM4, SPDBRK.C0M4) 
C0M_ALPHA(4) = A-COM4 
COM_BANK(4) = B-COM4 
COM-THRUST (4) = THR0T.C0M4 
COM_SPDBRK (4) = SPDBRK-COM4 
C 

STATUS = PPL$WAIT-AT_BARRIER(BARRIER-CMD) 

C 

101 CONTINUE 

END 
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Appendix B 

Communication and Synchronization Subroutine 

This appendix presents the communication and synchronization subroutine. This subroutine allows the TMS 
executive program to pass information in and out of the aircraft simulations by means of the shared common 
block variables. The barriers in this subroutine allow the executive program to maintain synchronization of 
the simulations. 

SUBROUTINE TMS 
C 

C OUTPUT FROM AIRCRAFT SIMULATION 
C 

1 (ALFDG, BNKCUR, BETDG , CXX, CXY, CXZ, CYX, CYY, 

1 CYZ , CZX, CZY, CZZ, PHIDG, THEDG , PSIDG, GAMDG, 

1 AZ, MACH, SX, SY, H, XD, YD, HD, EO, El, E2, E3, 

1 PDG, QDG, RDG, PWDG, QWDG, RWDG, DSB, XPCA, XPCS, 

1 PCR, T, TT, VT, 

C 

C INPUT FROM DECISION LOGIC 
C 

1 ALFCOM .BNKCOM, CSB, CPR) 

C EXTERNAL DEFINITIONS 
INTEGER* 4 PPL$GET_INDEX 

INTEGER*4 PPL$CREATE .BARRIER, PPL$WAIT_AT_BARRIER 
INTEGER*4 LIBSSTOP, LIB$PUT .OUTPUT 
C LOCAL DATA 

REAL MACH 

INTEGER*4 STATUS, MY .INDEX 
INTEGER*4 0NE.PAGE 

PARAMETER (ONE .PAGE = 512) 

C DATA FOR SHARING 

BYTE FRONT-GUARD (ONE .PAGE) 

INTEGER COPIES 

REAL AIRSPEED (4) , ALPHA (4) , BANKWND(4) , BETA (4) 

REAL DIRC0S(4,9) , EULER(4,3), GAMMA (4) 

REAL GL0AD(4) , MCH(4) , P0SITI0N(4,6) 

REAL QUAT(4,4) , R0TRATES(4,6) , SPDBRAKE(4) 

REAL STKRUD(4,3) , TIME(4), THRUST(4) 

C 

REAL C0M.ALPHA(4) , COM _BANK(4) , C0M.SPDBRK(4) , C0M_THRUST(4) 

BYTE REAR-GUARD (ONE .PAGE) 

C PUT SHARED DATA IN TO COMMON BLOCK 
COMMON /SHARED .DATA/ FRONT-GUARD, 

1 COPIES, 

1 AIRSPEED, ALPHA, BANKWND , BETA, 

1 DIRCOS, EULER, GAMMA, 

1 GLOAD, MCH, POSITION, 

1 QUAT, ROTRATES, SPDBRAKE, 

1 STKRUD, TIME, THRUST, 

C 

1 C0M.ALPHA, COM .BANK, 
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1 

1 


COM_SPDBRK , COM_THRUST, 
REAR_GUARD 


C 

STATUS = PPL$CREATEJ3ARRIER (BARRIER _DATA, ' BARRIER JDATA ’ , 
1 '/,REF (COPIES+1) ) 

STATUS = PPL$CREATE_BARRIER(BARRIER_CMD, i BARRIER_CMD 1 , 

1 %REF (COPIES+1) ) 

STATUS = PPL$GET_INDEX (MY_INDEX) 

C 

C****** PASS DATA TO TMS EXECUTIVE********* 

C 


AIRSPEED (MY_INDEX) 

= 

VT 

ALPHA (MY.INDEX) 

= 

ALFDG 

B ANKWND (MY .INDEX) 

= 

BNKCUR 

BETA (MY.INDEX) 

= 

BETDG 

DIRCOS(MY_INDEX, 1) 

= 

cxx 

DIRCOS (MY.INDEX, 2) 

= 

CXY 

DIRCOS (MY.INDEX, 3) 

= 

CXZ 

DIRCOS (MY.INDEX, 4) 

= 

CYX 

DIRCOS (MY.INDEX, 5) 

= 

CYY 

DIRCOS (MY.INDEX, 6) 

= 

CYZ 

DIRCOS (MY.INDEX, 7) 

= 

CZX 

DIRCOS (MY.INDEX, 8) 

= 

CZY 

DIRCOS (MY .INDEX, 9) 

= 

CZZ 

EULER (MY .INDEX , 1 ) 

= 

PHIDG 

EULER (MY .INDEX, 2) 

= 

THEDG 

EULER (MY .INDEX, 3) 

= 

PSIDG 

GAMMA (MY.INDEX) 

= 

GAMDG 

GLOAD (MY.INDEX) 

= 

AZ 

MCH (MY.INDEX) 

= 

MACH 

POSITION (MY.INDEX, 1) 

= 

SX 

POSITION (MY.INDEX, 2) 

= 

SY 

POSITION (MY.INDEX, 3) 

= 

-1.* H 

POSITION (MY.INDEX, 4) 

= 

XD 

POSITION (MY.INDEX, 5) 

= 

YD 

POSITION (MY.INDEX ,6) 

= 

-1. * HD 

QUAT(MY_INDEX, 1) 

= 

EO 

QUAT (MY.INDEX , 2) 

= 

El 

QUAT (MY.INDEX , 3) 

= 

E2 

QUAT (MY.INDEX, 4) 

= 

E3 

ROTRATES (MY.INDEX, 1) 

= 

PDG 

ROTRATES (MY.INDEX, 2) 

= 

QDG 

ROTRATES (MY.INDEX , 3) 

= 

RDG 

ROTRATES (MY.INDEX ,4) 

= 

PWDG 

ROTRATES (MY.INDEX , 5) 

= 

QWDG 

ROTRATES (MY.INDEX , 6) 

= 

RWDG 

SPDBRAKE ( MY _I NDEX ) 

= 

DSB / 60.0 

STKRUD (MY.INDEX, 1) 

= 

PCA 

STKRUD (MY.INDEX, 2) 

= 

PCS 

STKRUD (MY.INDEX, 3) 

= 

PCR 

TIME(MY. INDEX) 

= 

T 
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THRUST (MY_INDEX) = TT 

C 

STATUS = PPL$WAIT_AT_BARRIER (BARRIER .DATA) 
C 

STATUS - PPL$WAIT .AT J3ARRIER (BARRIER.CMD) 

C 

C* * * * * * * ACCEPT COMMANDS FROM EXECUTIVE******* 
C 

ALFCOM = COM.ALPHA (MY.INDEX) 

BNKCOM = COM.BANK (MY. INDEX) 

CSB = COM.SPDBRK (MY.INDEX) 

CPR = COM .THRUST (MY .INDEX) 

C 

RETURN 

C 

END 
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Appendix C 

Example Trim File 

This appendix presents an input or trim file for defining the initial conditions for a simulated aircraft. 


PHYSICAL 

CONSTANTS : 



31665.0 

WT 

(LBS) 

-WEIGHT 


22337.0 

IXX 

(SLUG*FT**2) -INERTIA ABOUT X 

AXIS 

120293.0 

IYY 

(SLUG*FT**2) -INERTIA ABOUT Y 

AXIS 

1 38945. C 

IZZ 

(SLUG*FT**2) -INERTIA ABOUT Z 

AXIS 

-2430.0 

IXZ 

(SLUG*FT**2) -XZ PLANE INERTIA PRODUCT 

457.3 

FSCG 

(IN) 

-FUSELAGE STATION CG 


0.0 

BLCG 

(IN) 

-BUTTOCK LINE CG 


101.6 

WLCG 

(IN) 

-WATER LINE CG 


687.5 

XNRF 

(IN) 

-X THRUST CENTERLINE 


18.9 

YNRF 

(IN) 

-Y THRUST CENTERLINE 


100.0 

ZNRF 

(IN) 

-Z THRUST CENTERLINE 



FLIGHT CONDITIONS AND MODELING OPTIONS: 
F LTHVEC » TRUE= THRUST VECTOR ON 
T LFCS » TRUE= FLT CONTROL SYTEM ON 
T LTHDMS » TRUE= DMS PLA SCHEDULE 
T LRTE » TRUE= R/T EQV AERO 
0.90000 MACHTR (N.D.) 

10000.00 HIC (FT) 

5000.00 X IC (FT) 

0.00 Y IC (FT) 

0.00 MUDGTR (DEG) 

1.00000 GLOAD (G) 

TCASE : 

1 

TRIM DRIVER VALUES: 

NXTR -NUMBER OF DRIVER VARIABLES 


4 


ELEMENT LIMITS 

VARIABLE 

NAME, UNITS 

3 0.000 

1.000 

ALFTR 

(RADIANS) 

7 -1.000 

1.000 

THETR 

(RADIANS) 

14 -2.500 

5.000 

PCSTR 

(INCHES ) 

16 31.000 

130.000 

DPSYTR C/. POWER) 

TRIM DRIVEN 

VALUES : 




NYTR -NUMBER OF DRIVEN VARIABLES 


4 

ELEMENT VARIABLE NAME, UNITS 
1 UD (FT/SEC2) 

3 WD (FT/SEC2) 

5 qD (RAD/SEC2) 

7 GAMZR (RADIANS) 

NFSY -OLD VAR RETAINED FOR FILE COMPATIBILITY 
0 

NFAS -OLD VAR RETAINED FOR FILE COMPATIBILITY 
0 

INITIAL CONDITIONS: 


24 



0.90000000 

MACHTR 

(N.D. ) 

O.OOOOOOOOE+OO 

BETTR 

(RADIANS) 

0.45296673E-01 

ALFTR 

(RADIANS) 

O.OOOOOOOOE+OO 

PIC 

(RAD/SEC) 

O.OOOOOOOOE+OO 

QIC 

(RAD/SEC) 

O.OOOOOOOOE+OO 

RIC 

(RAD/SEC) 

0 . 45298599E-01 

THETR 

(RADIANS) 

0 . OOOOOOOOE+OO 

PHITR 

(RADIANS) 

3 . 14159265E+00 

PSITR 

(RADIANS) 

0 . OOOOOOOOE+OO 

GAMTR 

(RADIANS) 

0 . OOOOOOOOE+OO 

DTVL 

(DEGREES) 

0 . OOOOOOOOE+OO 

DTVR 

(DEGREES) 

0.50865169E-05 

PCATR 

(INCHES ) 

0 . 24973108E-01 

PCSTR 

(INCHES ) 

0 . 61914313E-03 

PCRTR 

(LBS ) 

77 . 252785 

DPSYTR 

('/. POWER) 

O.OOOOOOOOE+OO 

DPASTR 

(*/, POWER) 

-0 . 18907314E-01 

DSSYTR 

(DEGREES) 

-0 . 12405217E-05 

DSASTR 

(DEGREES) 

O.OOOOOOOOE+OO 

DASYTR 

(DEGREES) 

-0.31018224E-05 

UD 

(FT/SEC2) 

O.OOOOOOOOE+OO 

DRSYTR 

(DEGREES) 

-0.15484100E-04 

DRASTR 

(DEGREES) 

3.4465680 

DNSYTR 

(DEGREES) 

0 . OOOOOOOOE+OO 

DNASTR 

(DEGREES) 

3.6334312 

DFSYTR 

(DEGREES) 

-0 . 95367432E-06 

DFASTR 

(DEGREES) 

O.OOOOOOOOE+OO 

CSB 

(DEGREES) 


$END OF DATA READING SECTION 
CASE SELECTIONS: 

TCASE 1 STRAIGHT & LEVEL STEADY STATE 

2 COORDINATED TURN STEADY STATE 

3 PULL-UP STEADY STATE 
TRIM VALUE SELECTIONS: 


TRIM 

DRIVER 

ARRAY 


TRIM OUTPUT ARRAY 

1 

MACHTR 

(N.D.) 

1 

UD 

(FT/SEC2) 

2 

BETTR 

(RADIANS) 

2 

VD 

(FT/SEC2) 

3 

ALFTR 

(RADIANS) 

3 

WD 

(FT/SEC2) 

4 

PIC 

(RAD/SEC) 

4 

PD 

(RAD/SEC2) 

5 

QIC 

(RAD/SEC) 

5 

QD 

(RAD/SEC2) 

6 

RIC 

(RAD/SEC) 

6 

RD 

(RAD/SEC2) 

7 

THETR 

(RADIANS) 

7 

GAMZR 

(RADIANS) 

8 

PHITR 

(RADIANS) 

8 

PHIZR 

(RADIANS) 

9 

PSITR 

(RADIANS) 

9 

THE 

(RADIANS) 

10 

GAMTR 

(RADIANS) 

10 LAMDA 

(RADIANS) 

11 

DTVL 

(DEGREES) 

11 

fytot 

(G S) 

12 

DTVR 

(DEGREES) 




13 

PCATR 

(INCHES ) 




14 

PCSTR 

(INCHES ) 




15 

PCRTR 

(LBS ) 




16 

DPSYTR 

(*/. POWER) 






17 

DPASTR 

(*/. POWER) 

18 

DSSYTR 

(DEGREES) 

19 

DSASTR 

(DEGREES) 

20 

DASYTR 

(DEGREES) 

21 

DAASTR 

(DEGREES) 

22 

DRSYTR 

(DEGREES) 

23 

DRASTR 

(DEGREES) 

24 

DNSYTR 

(DEGREES) 

25 

DNASTR 

(DEGREES) 

26 

DFSYTR 

(DEGREES) 

27 

DFASTR 

(DEGREES) 

28 

CSB 

(DEGREES) 
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Table I. Summary of Weight, Center of Gravity, and Inertia 


Weight, 

lb 

Center-of-gn 

ivity locations 

Moments and product of inertia, slugs/ft 2 | 

Fuselage 

station, 

in. 

Water 

line, 

in. 

J XX 

I Y y 

hz 

Ixz 


TV aircraft 


33310 

455.0 

102.8 

23000 

151293 ! 

169 945 

-2971 


1 

Baseline aircraft 


| 31 665 

457.3 j 

101.6 

22 337 

120293 

138 945 1 

-2430 
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Table II. Dimensional Data 0 


Total airplane: 

Net wetted area (minus engine nozzles) , ft 2 

Overall length, ft 

Overall height, ft 

Wing: 

Area, Sref> ft 2 

Wetted area (including launchers and aileron actuator fairings), ft 

Span, 6 ref> ft 

Aspect ratio 

cref. ft 

Leading-edge sweep, deg 

c/4 sweep, deg 

Taper ratio 

Dihedral, deg 

Leading-edge flaps: 

Deflection (positive leading edge down), deg — 

Maneuvering 

Takeoff and landing 

Differential 

Actuator 18 deg/sec rate limit 

Trailing-edge flaps: 

Deflections (positive trailing edge down), deg— 

Takeoff 

Landing 

Actuator 18 deg/sec rate limit 

“From reference 13. 


2028 

56.0 

15.3 


400 

562 

37.42 

3.5 

11.52 

26.7 

20 

0.35 

-3 


. . 0, 34 

. . 12, 34 
. . . . ±3 

l/(s/20 + 1) 


. +17, +30 

. +17, +45 

l/(s/20 + 1) 
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Table II. Concluded 


Ailerons: 

Deflections (positive trailing edge down), deg — 

Takeoff and landing 

Maneuvering 

Actuator 100 deg/sec rate limit 

Horizontal tails (HT): 

Exposed area, ft 2 

Aspect ratio 

c/4 sweep, deg 

Span, EF> ft 

QiTi ft 

Deflections (positive trailing edge down), deg — 

Symmetric 

Maximum 

Actuator 40 deg/sec rate limit 

Vertical tails (VT): 

Area, ft 2 

Wetted area, ft 2 

c/4 sweep, deg 

Cant (tip out), deg 

^VT> ft 

Tail length (c/4 to cy x /4), ft 

Rudders: 

Deflection, deg 

Actuator 61 deg/sec rate limit 

Speedbrake: 

Planform area, ft 2 

Span, ft 

Chord, ft 

Maximum deflection, deg 


. -25, +45 
. -25, +25 
l/(j/48+l) 

. . . 88.1 
... 2.4 

. . . 42.8 
. . 14.67 

. . . 6.28 


. . -24, +8 
. -24, +10.5 
l/(s/30 + 1) 


52.0 each 
104.0 each 
. . . 35.0 
. . . 20 
. . . 6.99 
. . 10.18 


. ... ±30 
l/(s/40 + 1) 

. . . 13.9 
... 2.5 

. . . 5.57 
... 60 
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Tabic III. Time Required by TA To Perform o Captures 
[All runs started at h = 25000 ft and had ±2° capture criteria] 


Initial 

Final 

Initial 

Average time 

Minimum time 

Time by 

Maximum 

a, deg 

a, deg 

M 

by pilot, sec 

by pilot, sec 

TA, sec 

overshoot, deg 




Baseline 

aircraft 



4.4 

30.0 

0.60 

5.12 

4.35 

1.91 

1.9 

4.4 

40.0 

.60 

2.88 

2.30 

“2.28 

2.4 

23.5 

30.0 

.30 

4.93 

3.78 

1.00 

1.4 

23.5 

40.0 

.30 

6.56 

5.95 

1.81 

1.6 

10.0 

0.0 

.40 

2.50 

1.99 

1.34 

1.0 

20.0 

0.0 

.32 

5.86 

5.25 

1.88 

2.0 

30.0 

0.0 

.27 

7.06 

5.68 

2.38 

2.0 

TV aircraft 

4.4 

30.0 

0.60 

4.70 

3.84 

1.09 

1.7 

4.4 

40.0 

.60 

4.45 

3.46 

2.97 

2.6 

4.4 

50.0 

.60 

4.76 

5.31 

2.41 

.2 

23.5 

30.0 

.30 

2.11 

1.09 

.81 

1.2 

23.5 

40.0 

.30 

2.69 

1.41 

1.38 

1.2 

23.5 

50.0 

.30 

3.39 

1.79 

1.78 

1.6 

10.0 

0.0 

.40 

2.18 

2.18 

1.12 

.4 

20.0 

0.0 

.32 

2.11 

1.66 

1.60 

.7 

30.0 

0.0 

.27 

4.60 

4.54 

1.89 

.6 


“Capture criteria relaxed to ±2.4°. 


Table IV. Time Required by TA To Perform 90° p Captures 
[All runs started at h = 25 000 ft, Initial p = 0°, and Final p = 90°] 


Initial 

Capture 

Average time 

Minimum time 

Time by 

Maximum 

a, deg 

criteria, deg 

by pilot, sec 

by pilot, sec 

TA, sec 

overshoot, deg 

Baseline aircraft 

10 

±5 

4.10 

3.07 

1.43 

3.8 

20 

±8 

8.90 

6.70 

4.90 

6.0 

TV aircraft 

10 

±5 

2.15 

1 1.47 

1.53 

2.8 

20 

±5 

5.00 

4.40 

2.22 

2.7 

30 

±5 

5.17 

2.75 

2.50 

3.9 
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Figure 1. Configuration of aerodynamic surfaces, definitions of axes, and sign convention (ref. 14). 


Vane size 



Figure 2. TV system. Vane cant angle 48°; maximum vane deflection ±30°; maximum deflection rate 
80 deg/sec. (All linear dimensions are in inches.) 
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g limiter 


Position 

limiter 



Figure 3. P = Signal. Longitudinal control augmentation system (ref. 20). 


















Figure 4. Lateral and directional control augmentation system (ref. 20). 
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Figure 5. TDG TA aircraft system. 



(a) Longitudinal command system. 



(b) Lateral command system. 
Figure 6. Block diagram of TA system. 
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- 40.0 - 30.0 - 20.0 



(c) Angle-of-attack and stabilator responses up to (d) Angle-of-sideslip and rudder responses up to 
a =-10°. (3= 20°. 

Figure 7. Departure of baseline aircraft with full lateral stick input. Initial M — 0.30; h — 10000 ft. 
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Rudder 









.00 3.00 4 00 5.00 T 0 1.00 2.00 3.00 4.00 

Time, sec Time, sec 


(a) Initial M = 0.60, baseline aircraft. 


o 



Time, sec Time, sec 


(b) Initial M = 0.30, baseline aircraft. 

Figure 8. Performance of TA for longitudinal captures. 







a, deg a, deg 

- 1 0.0 10.0 30,0 50.0 70.0 90.0 ^ - 10.0 10.0 30.0 50.0 70.0 90.0 




Time, sec 


Initial M = 0.40, baseline aircraft. 


(d) Initial M = 0.32, baseline aircraft. 



8 > 

T3 

6 


o 



Time, sec 


(e) Initial M = 0.27, baseline aircraft. 


(f) Initial M = 0.60, TV aircraft. 


Figure 8. Continued. 






Actual a 
Commanded a 


.00 3.00 4.00 5.00 

Time, sec 



(g) Initial M = 0.60, TV aircraft. 



(h) Initial M = 0.30, TV aircraft 
Figure 8. Continued. 





006 ooz oos ooe ooi om- 006 ooz oos ooe ooi o - oi- 

6ep 6ep *x> 


o 




i) Initial M = 0.30, TV aircraft. 


(j) Initial M = 0.40, TV aircraft. 



o 



Time, sec 


(k) Initial M = 0.32, TV aircraft. 


(1) Initial M = 0.27, TV aircraft. 


Figure 8. Concluded. 



















(e) a = 30 


Figure 9. 
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Aircraft simulation 
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DATA 
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Com and sync ^ 

r Ti 

BARRIER 



\( 11 
BARRIER. 
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DATA 
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CMD 
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Figure 12. Operational flow between TMS executive program and one aircraft simulation. Additional aircraft 
are simulated by concurrent execution of multiple copies of aircraft simulation model. 
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Figure 14. 


Simulation of multiple aircraft; positions plotted at 1-sec intervals. 










o 



(a) a response. 



Time, sec 


(b) Wind-axis bank-angle response. 

Figure 16. Commanded and actual responses of guided aircraft. 
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